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procedure COMPUTEMlNERROR( Aggregate x, Aggregate^, integer 1) 

1 . if subTree[x, y, I) .computed = true 

2. return [subTree[;c, y, /]. error, subTree[x, y, /]. aggregates] 

3. niinError := minErrorl := minError2 := co 

4. if a* is a leaf { 

5. minErrorl := 2 se sZ)0, /) * (lsp(s, x, {y}, W A ) - lsp(s, x)) 

6. if/>0 

7. minError2 := E se s£>0, /) * (lsp(s, x, {*}, W A ) - lsp(s, x)) 

8. if minErrorl < minError2 

9. [subTree[x, y s /]. error, subTree[x, y, /]. aggregates] := [minErrorl, 0] 

10. else 

11. [subTree[x, y, /]. error, subTree[x, y, /]. aggregates] := [minError2, {x}] 

12. } 

13. if x has a single child u { 

14. [minErrorl, aggregates 1] := COMPUTEMINERROR(m, y, I) 

15. if/>0 

16. [minError2, aggregates2] := COMPUTEMlNERROR(w, x, I - 1) 

17. if minErrorl < minError2 

18. [subTreefx, y /].error, subTree[x, y, /]. aggregates] := [minErrorl, aggregatesl] 

19. else 

20. [subTree[x, y, /]. error, subTree[x, y t /] .aggregates] := [minError2, aggregates2 u {x}] 

21. } 

22. if x has children u and v { 

23. for/ ~0to/ { 

24. [minErrorl, aggregatesl] := COMPUTEMlNERROR(w, y, i) 

25. [minError2, aggregates2] := C0MPUTEMinErR0R(v, y, k - i) 

26. if minErrorl + minError2 < minError 

27. minError := minErrorl + minError2 

28. aggregates := aggregatesl kj aggregates2 

29. } 

30. for/:=0to/-l { 

31. [minErrorl, aggregatesl] := COMPUTEMlNERROR(w, x, i) 

32. [minError2, aggregates2] := COMPUTEMlNERROR(v, x, k-i-l) 

33. if minErrorl + minError2 < minError 

34. minError := minErrorl + minError2 

35. aggregates := aggregatesl u aggregates2 ^ {x} 

36. } 

37. [subTree[x, y, /]. error, subTree[x, y f /]. aggregates] := [minError, aggregates] 

38. } 

39. subTree[x, y t /]. computed := true 

40. return [subTree[x, y t /]. error, subTreefx, y, /]. aggregates] 
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procedure COMBINEMlNERROR() 

1 . for i : = 1 to m 

2. for j = 0 to k { 

3. Tit/']. [error, aggregates] := COMPUTEMlNERROR(r(7i ), e,y) 

4. X\\j\. [error, aggregates] := [ oo 5 0] 

5. } 

6. for j = 0 to k 

7. X ,[/']. [error, aggregates] := Ti [/'].[ error, aggregates] 

8. for i = 1 to m 

9. for; =OtoA: 

10. for / = 0 to j 

11. if (JTj.i [/] .error + T,[j - /] .error < Xi[f] .error) { 

1 2 . Xi\f[ . error = JQ. i [/] . error + T t \J - /] . error 

13. Xi\j]. aggregates = X iA [/]. aggregates u W - /].aggregates 

14. } 
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procedure COMPUTEWEIGHTSCUMULATIVE () 

1. for each b e B, set W min (b) : = 0 

2. for / :=1 to r { 

3. W:=W min 

4. Choose a random subset R cz 5, of ABRs 

5. for each b e set ^(6) to a random weight in [0, L] 

if Hsese(s 9 B(s, W)) < Z se se (s, B(s, W min )) 
W min := W 



} 

return PT« 



i^TC. 9 

procedure ComputeWeightsMax(Q) 
1 .for each b e .B, set Wold(b) := 0 
2. while (P£ 2 B 

1 Wold(b) < ( 

jS/j*G£ij-l) 

2 ) *lspmax) f3. Let 

QO be a new set of inequalities that result when the value Wold(b) is 
substituted for each variable W (b)only on the LHS of each inequality in 
Q 4. Set Wnew(b) to the smallest possible value such that each 
inequality in QO is satisfied when Wnew(b) is substituted for variable W 
(b) in QO 5. if Wnew= Wold 6. return Wnew 7. else 8. Wold := Wnew 
9.g 10. return "there does not exist a weight assignment W n 
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procedure COMPUTEWElGHTSTwoABR() 

1 . Set V opt := vOO, E := E opt := Hs^se(s y b 0 

2. for j:=l to n { 

3. E := E + e(sj 9 62) - e(sj> b\) 

4. \iE<E opt 

5. V opt := v(s' } +\) 9 E 0 pt = E 

6. } 

7. return V ODt 



